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Nomenclature 


Note: Units are used in which c = h = 1. 
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speed of light 

Planck constant divided by 2ir 

as a subscript, refers to lab frame 

as a subscript, refers to center of mass frame 

flux factor for masses i and j 

energy of particle i 

3- momentum of particle i 

4- momentum of particle i 
4-momentum of virtual particle j 
invariant amplitude 

coupling constant 

coupling constant between particles 1, 3 and exchange particle x 

pion-nucleon coupling constant 

pion-nucleon-delta coupling constant 

electromagnetic coupling constant 

mass of particle j 

mass of the A particle 

mass of the pion 

mass of the exchange particle (e.g. x = n, the pion) 

mass of the proton 

mass of a nucleon 

total decay width or rate 

partial decay width for exclusive decay 4 — > 56 
differential decay width 
statistical factor 
total cross section 

total cross section due to direct term 
total cross section due to exchange term 
total cross section due to interference term 
infinitesimal element of differential cross section 
invariant differential cross section 
angular distribution 

spectral distribution in terms of energy of particle i 
spectral distribution in terms of kinetic energy of particle j 
a nucleon 

lab angle of particle 3 with respect to particle 1 
cm angle of particle 4 with respect to particle 1 
maximum angle of particle 3 with respect to particle 1 
mandelstam variables 
value of variable t at 6 = 0 
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d&2 

d® n 

dLips 


value of variable t at 9 = n 
value of variable u at 0 = 0 
value of variable u at 6 = vr 

= ^threshold) value of variable s at the reaction threshold 

the neutral pion 

the A baryon of charge +1 

the photon 

mass energy unit, one billion electron volts 

millibarns 

picobarns 

steradian, unit of solid angle 

phase space factor for 2 particle final state 

phase space factor for n particle final state 

Lorentz invariant phase space factor 

mean lifetime 

one pion exchange 


Em 2 = Yli = i m i = m i + m 2 + + m| 
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Abstract 


A one pion exchange scalar model is used to calculate differential and total cross sec- 
tions for pion production through nucleon-nucleon collisions. The collisions involve 
intermediate delta particle production and decay to nucleons and a pion. The model 
provides the basic theoretical framework for scalar field theory and can be applied to 
particle production processes where the effects of spin can be neglected. 

1 Introduction 

An effort is underway to include meson production in the space radiation transport code known 
as HZETRN. Previous work used high energy meson production models based on scaling [1, 
2, 3, 4, 5]. These models yield nucleon, pion, and other particle production cross sections. 
Basing the transport codes instead on fundamental models of hadronic interactions holds out 
the prospect of providing unified and more versatile descriptions of the interactions of interest. 
The quark - gluon model is the most fundamental, but its complexities are such that calculations 
are more of a qualitative kind [6]. The present work is based on the meson exchange model of 
hadronic (strong) interactions [6]. The model may be considered an effective quantum field 
theory in which the exchange particle is the pion. In the model’s simplest form, all particles, 
nucleons and exchange particles alike, are treated as scalar, spin zero particles. This version 
of the model makes no distinction between proton and neutron other than by mass, but since 
the masses are nearly equal, the mass of the proton is used throughout as the nucleon mass. 
As a further simplification, only the pion is used as the exchange particle. This scalar one 
pion exchange (OPE) model is used to determine the amplitudes, to first and second order in 
coupling constants, of A decays and nucleon - nucleon ( NN ) interactions producing 2 - body 
and 3 - body final states. The amplitudes are inserted into standard formulas derived from 
Fermi’s golden rule to give differential decay widths for A — + Nir and differential cross sections 
for the elastic NN — > NN and inelastic NN —> IV A interactions. Manipulation of the phase 
space factor in the differential cross section leads to an expression for the invariant distribution, 
a differential cross section that is a function of Mandelstam variable t. Changes of variable 
lead to angular distributions in the center of momentum (cm) frame and angular and spectral 
distributions in the lab frame. Of the various distributions, the invariant distribution reveals 
itself to have a particularly simple form. Allowed ranges of variables are found, and integration 
over the variables yields total cross sections. 

2 Scattering for 2 - body final states 

In 2 - body decay, the lowest order term in the scattering matrix has one vertex, and thus is 
first order in the coupling constant g n NN ■ In 2 - body scattering, the lowest order term in the 
scattering matrix has two vertices, and thus is second order in the coupling constant. 2 - body 
reactions will be represented as 1 + 2 — ► 3 + 4, where the numbers represent particles. Both 
decay and scattering processes have a 2 - body final state. In this section, the basic formalism 
for calculating 2 - body final state cross sections is presented. 
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2.1 Scattering amplitude from Feynman rules 

In this subsection, it will be shown how to generate a scattering amplitude from a set of Feynman 
rules. This amplitude is the basic input to physical observables, such as cross sections and decay 
rates. The Feynman rules used below apply to a set of particles that are scalars. Griffiths [7] 
labels the particles as A,B,C each with different masses. The exchange particle is particle C. 
We consider the Feynman rules for the scalar ABC theory of Griffiths [7]. The rules of Peskin 
and Schroeder [8] and Griffiths [7] are the same except that Griffiths has rules leading to the 
invariant amplitude — iM , whereas Peskin and Schroeder have the rules giving +iM. This 
discrepancy is due to different definitions relating the S-matrix to the invariant amplitude. We 
follow the convention of Peskin and Schroeder in the present work. Another difference is that 
Peskin and Schroeder always include a factor ie in the denominator of propagators. We include 
this factor in the form imT in the case of A decay to account for the decay width. The Feynman 
rules are now discussed. The quantity 


iM{2-jr) 4 5 4 (p 1 +p 2 ~P3 ~P4 Pn) , (1) 

is given by the following rules [7, 8] for scalar ABC theory [7]. 

1. Notation. Label the incoming 4-monrenta as p\ and p 2 , and the outgoing momenta as 
P3 , P4 ' ' ' Pn ■ Label the internal 4-monrenta as q \ , q 2 . . . . Put arrows on each line to conserve 
current, which identify incoming and outgoing particles. 

2. Coupling constant. For each vertex, write down a factor 

-ig 

where g is the coupling constant. (In our simple scalar theory, this has dimensions of GeV, 
but in realistic theories g is always dimensionless. In the theory of Quantum Electrody- 
namics, g is the fine structure constant a = 1/137.) 

3. Propagator. For each internal line, write a factor 

i 

qj ~ m J + 

where qj is the 4-momentum of the internal particle and nrij is its mass. Note that qj =/=■ rrij , 
because the internal particle is virtual and does not obey the Einstein relation, p 2 = m 2 . 
We say that such a particle is off mass shell. Also, as mentioned above, the factor ie is 
omitted except in the case where decay width is significant. 

4. Conservation of 4-momentum. At each vertex, write a delta function of the form 

(27r) 4 5 4 (fci + k 2 + k 3 ) 

where the k's are the 4-monrenta flowing into the vertex. If the particle flows out, then k 
should have a negative sign in front. 
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5. Integrate over internal momenta. For each internal line, write down a factor 


1 



and integrate over all internal momenta. 


2.2 Decay width from Feynman rules 

Most baryons decay into a final set of particles. Therefore, it is important to be able to describe 
this decay. This subsection shows how to calculate the decay width from Feynman rules. The 
Feynman rules may be used to calculate the decay width of particle 3 decaying, as in the reaction 

3— >1 + 2. (2) 

The Feynman diagram for this reaction is shown in figure 1. This diagram is particularly easy 
to evaluate since there are no internal lines. Applying the Feynman rules to the lowest order 
diagram (see Griffiths [7]), gives 

iA/f(27r) 4 (S 4 (pi -p 2 -p 3 ) = (-ig)(2ir) 4 5 4 (p 1 -p 2 -p3), (3) 


and canceling terms gives 


M = —g . 


(4) 



Figure 1: Lowest order contribution to the decay 3 — > 1 + 2. 


Step 5 in section 2.1 is not needed because there are no internal particles for the first order 
decay. Substituting equation (4) into equation 6.30 of Reference [7] gives 




t 167rmo 


\J m\ + + wij — 2m\ — 2 m\rn^ — 2 


(5) 


where g = g nN A- 
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2.3 Scattering via massive exchange particle 

In this section, we consider scattering by means of an exchange particle with mass m > 0. 
Since the exchange particle is a virtual particle, its mass is allowed to be off-shell. In the case 
of the pion, this means that the square of the 4-momentum transferred by the virtual pion 
t = (pi — P3) 2 is in general not equal to the squared mass of a real pion, t / m 2 . In treating 
the reaction NN — > IV A, the ABC theory is simply expanded to include two types of vertices, 
along with their respective coupling constants: NNtt and NAir. 

2.3.1 Amplitude 

Consider the reaction 


1 + 2 — > 3 + 4 , 


with the diagram shown in figure 2. The particles have masses mi, m 2 , m 3 and m 4 , respectively, 
and the mass of the intermediate exchange particle is rn x . Applying the Feynman rules to the 
lowest order diagram (see Reference [7]), gives the direct amplitude as 


iMd(27r) 4 <f(pi +P2-P3 - pa) 


— (— i) 2 9l3x92Ax 

= (~ ■ i) 2 9l3x92Ax 


d 4 q i 
(27 r) 4 q 2 - 777.2 


(27r) 4 <P(pi -P 3 ~ <?)(27 t) 4 <5 4 (<7 + P2- Pa) 


('P\ -P3) 2 ~ m 2 


(27 t) 4 <P(pi -P3 + P2- Pa) , 


and canceling terms gives 


where t = {p\ — P3) 2, ■ 


M d = 


9l3x92Ax 
t — m 2 


( 6 ) 

( 7 ) 



Figure 2: Lowest order contribution to the reaction 1 + 2 — > 3 + 4. 


The exchange term is shown in figure 3 and is given by 


M e 


gl4^g23x _ ffl4rff23a: 
( pi — pa) 2 — m 2 u — ml 


( 8 ) 
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where u = (pi — Pa) 2 ■ The plus sign of the exchange term in equation (8) is consistent with the 
scalar theory in which all particles are spin 0 bosons with wave functions that commute. Had 
spin been considered, then for the case of nucleon scattering, in which the nucleons are fermions 
with spin 1/2, the exchange term in equation (8) would have taken on a minus sign reflecting 
the fact that the wave functions of the two fermions anti-commute [8] (p. 119). 



Figure 3: Exchange term of reaction 1 + 2 — ► 3 + 4. 

For the reaction NN — ► NN, the direct and exchange terms are physically indistinguish- 
able in so much as the energies and exit angles of a final state particle are the same for both 
terms. According to the rules of quantum mechanics, the probability (i.e. cross section) to 
produce the final state particles is determined by first adding the amplitudes of the two terms, 
and then squaring the complete amplitude, (rather than the reverse, as would be the case for 
distinguishable final states [9] (p. 10). 


Table 1: Values of t and u for 0 3 i a b = 0.4 


E 3 iab(GeV) 

t(GeV 2 ) 

u(GeV 2 ) 

2.98 

- 2.07 

- 3.19 

1.30 

- 5.23 

- 0.035 


For each energy, there is a unique value of t that determines that energy. This value of t 
represents the direct term. There is also a unique value of u that determines that energy. See 
Table 1. This value of u represents the exchange term. The invariant amplitude AA consists of 
the sum of direct (t-channel) and exchange (u-channel) terms and is finally 


AA — Aid + AA e — 


ff!3xff24x 9l4x923x 
t — m% u — m%. 


(9) 
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The Mandelstam variables s,t,u satisfy 


4 

s + t + u = m\ + m\ + + m\ = ^ mf = Em 2 , 

i — 1 


where .s = (p± + P2) 2 ■ 

2.3.2 Invariant t distribution 

Final formulas for the Lorentz invariant angular cross sections will now be presented. 
[10] gives 


da S 
dt 167rAi2 


\m\ 2 = 


s 


167rAi2 


9l3x92Ax 9lAx923x 
t — mi u — ml 


with 


\ / 2 2\2 A 2 2 

A ij — (s — m, — TBj ) — rrij . 

The direct, exchange, and interference terms are 

d&d _ $ | ja |2 _ $ 9l3x92Ax 

dt 167rAi2 d 167rAi2 (t — mi) 2 ’ 


dc T e 
dt 
dc n 
dt 


S 


167rAi2 

cS 


\M e \ 2 = 

imi 2 = 


cS 


9\Ax923x 

167rAi2 (u — mi) 2 

2 - ^ 2 M d M e = 


•S 2 9l3x924x9lAx923x 


167rAi2 r M 167rAi2 " "" L 167rAi2 (t — m, 2 )(u — m 2 ) 
For future purposes, we write these more simply as 


da d 

I< 

dt 

(' t~m 2 ) 2 ’ 

da e 

K 

dt 

(u — mi) 2 

d<7 t 

2 K 

dt 

(t-m 2 )(u- 


’ 


with the definition 

K = K(s) = 


S 


167rAl2 


2 2 _ 
9l3x92Ax ~ 


$ 9l3x92Ax 


167t[(s — — m^) 2 — 4 mfm|] ’ 


which assumes 


9 23x — 9l3x , 
92Ax — QlAx • 
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( 10 ) 

Reference 

( 11 ) 

(12) 

(13) 

(14) 

(15) 

(16) 

(17) 

(18) 

(19) 

( 20 ) 
(21) 



This equivalence of coupling constants between the direct and exchange terms assumes particles 
1 and 2 are the same. This is the case for the following discussions in which the particles are 
both nucleons. K(s) has a singularity when (s — m\ — 77T. 2 ) 2 — 4 = 0, and it will turn out 
that this gets canceled by other terms. See the note following equation (111). 

Elastic scattering (e.g. NN — > NN) 


Consider the example of the reaction NN — > NN where the exchange particle is a pion, 


which gives 


9ld>x ~ 

924x z 

= 9NN-K — 9nNN , 

(22) 

777-1 = 

777-2 = 

■ 7773 = 7714 = 771JV j 

(23) 

m x = 

777-7T = 

: 777- , 

(24) 

s = 

1/2, 


(25) 

Al2 

= s(s 

- 4m^) , 

(26) 


and 


K = 


S 


167rAi2 


2 2 

~9l3x924x ~ 


9tcNN 


327 ts(s — 4 m 2 N ) 


(27) 


Inelastic scattering (e.g. NN — > 1VA) 


Consider the example of the reaction NN -» 1VA where the exchange particle is a pion, 


which gives 


9l3x 

924x 

mi 

777-4 

m. x 

s 


9NNit = 9-kNN , 
9NAn = 9irNA , 
777-2 = = 777-AT , 

777-A , 

777-Tr = 777 , 

1 , 


A 12 = s(s - 4m jy) , 


(same as elastic case), 


and 


K = 


S 


167rAi2 


2 2 _ 
~9l3x924x ~ 


2 2 
9-KNN9-KNA 

167t.s(s — 4m^) 


(28) 

(29) 

(30) 

(31) 

(32) 

(33) 


(34) 


(35) 
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2.3.3 Spectral distribution in the lab frame 

The spectral distributions da/dE^i and da/dT^ are obtained from the invariant distribution 
(11) by expressing t and u in terms of E 41 , the energy of particle 4 in the lab frame, or in terms 
of T 41 , the kinetic energy. Expanding 


t = (P2 -P4) 2 , 

and noting that in the lab frame E 2 = m 2 and P 2 = 0, gives 

t = m| + m\ — 2 m 2 E±i 


and 


u = Em 2 — s — m 2 — m 2 + 2 m 2 -E 4 / , 


with 


2 2 1 2 1 2 1 2 

Em = m 1 + m 2 + m 3 + ?n 4 . 

Both t and u are expressed in terms of T 41 by the substitution E 41 = T 41 + m 4 , as in 

t = m 2 + m\ - 2 m 2 (T 4 ; + m 4 ) , 
u = 'Em 2 — s — m 2 — m 2 + 2?B2(T 4 / + m 4 ) . 


(36) 

(37) 

(38) 

(39) 

(40) 

(41) 


The range of possible values of E 41 and T 41 are bound by their minimum and maximum values, 
which are in turn found by inverting the preceding formulas [10] and substituting to and t n for t 
from equation (56), which is listed later in this paper. The above formulas for t and u can now 
be used to generate da/dE 4 i and da/dT 4 i from da/dt and dcr/du. 


2.3.4 Angular distribution in the cm frame 

The angular distribution is formed by substituting the expression for At given by equation (9) 
[10] and by expressing t and u in terms of angle 0. For a 2 - body final state, this is done 
by writing the momenta P 3 and p 4 in terms of the initial energy in the form of s and A. The 
denominator in the direct term in equation (7) contains 


t = ( P4~P2f 


ml + ml - 2\J-pl + ml\jY >2 + m 2 + 2 |P4| |P 2 1 cos 9 

\/ p /cm + m 4 yj Pi cm + m 2 + 2 1 Pf cm \ \ Pi cm \ cos 6 (cm frame) 


m| + m| — 2 


ml + m\ + — 


- \J (A 12 + 4 .sm|)(A 3 4 + 4 smf) + \j A 12 A 34 cos 6 


(42) 



where p cm = pi = — P2 and p/ C m = P3 = — P4 are the initial and final 3-nromenta of the 
particles in the cm frame. The denominator in the exchange term in equation (8) contains 

u = (p 3 ~ p 2 ) 2 

= m\ + m\ — 2E 3 E 2 + 2p3.p2 

= + m 2 — 2E 3 E 2 — 2p4.p2 (because P3 + P4 = 0 in cm frame) 

= m§ + — + P2 + m 2 ~ 2 1 P4 1 1 P2 1 cos 6 (in cm frame) 

= m§ + ml -2\Jp} cm + ml\Jpicm + m 2 ~ 2 |P/em||Picm| cos 8 (in cm frame) 

= Em 2 — s — t , (43) 


where 8 is the angle between p2 and P4 in the cm frame. The angular distribution becomes [10] 

da 


dEl^ 


64lT 2 S V A12 

S I A34 




64ir' 2 s V A12 


9l2x9Mx 9l2x93Ax 
t — m 2 u — m 2 


with the definitions 


t = mi + mi + 


2s 


- \J (A12 + 4sm|)(A 34 + 4 sm\) + yj Ai 2 A 34 cos 8 


(44) 


(45) 


and 


u = Em 2 — s — t . (46) 

In the cm frame, the angular distributions of particles 3 and 4 are equivalent [10]. 

da/dQ 3c = da/d^l^c . (47) 

Comparing equation (11) with equations (44) and (45), the differential cross sections take a 
simpler form when written in terms of the Mandelstam variables s, t and u, instead of the 
momentum |p| and angle 8. 


2.3.5 Angular distribution in the lab frame 

The angular distribution in the lab frame may be expressed in terms of the angular distribution 
in the cm frame [10], 


da 

dQ 3 i 


4 m 2 s |PidlP3f| da 

VA12A34 | E u + m 2 - ^E 3i cos 8 i 3 i\ d ^i3c 


4:1712 S 
'J Al2 A34 


2/^ 

\E3i(rrf 4- + m\ 


- (El - m 2 3 ) 3 / 2 da 

m 2 + 2E 1 m, 2 ) - 2m\(Ei + m 2 )| dQ i 3c 


(48) 
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For computational purposes, both E 3 i and da/d£l\ 3c , which contain t, may be written as func- 
tions of $i 3 1 and the energy of particle 1 in the lab frame Eh (instead of s), by using the following 
relations for s(Eu) and t(En, 6 3 i), 

s(Eu) = m\ + m\ + 2Eim 2 , (49) 


and 


t(Ei,0 3 i) = ml - m\ + 2m 2 [E 3 (Ei,6 3 i) - Eu\ . 


(50) 


The expression for s is found by expanding s = (p\i + P21) 2 and noting that in the lab frame 
P 2 = 0 and E 21 = m 2 - The expression for t is found, following reference [10], by expanding 
t = (p2 — Pi ) 2 and eliminating E 41 using conservation of energy, Eh = En +m 2 — E 3 i. Then [10] 


ab ± J Ef, - m\ cos{0 3i ) J a 2 - m|[6 2 - ~[Ef t - m^)cos 2 (9 3 i)] 

Esi(Eu, e 3l ) = y— ,2 ,1 Z-2E - oTF ^ . (51) 


b 2 - (E'h - mf)cos 2 (9 3 i) 


where 


a = 


s + rn'l — m\ 


(52) 


and 


b = Eu + rn 2 . 


(53) 


Note the double valued nature of the energy expression (51). A discussion of how to handle this 
in computations is discussed in reference [10]. 


2.3.6 Singularities 

The differential cross sections in equations (16), (17), and (18) have singularities at 

t = m 2 , 


and 


u = m 2 , 


(54) 

(55) 


where m = m x . In this section, we show that these singularities are never encountered for 
physical values of t and u. The formulas derived in reference [10] for the extreme values of t 
denoted as to = t(9 = 0) and t n = t {6 = 7 r), are 


to{U) = ^ (ml - ml - m\ + m|) 2 - (y / Ai2 =F V^a) 2 


(56) 


The symbol T means that a negative sign is used for to and a positive sign is used for t n . A 
useful result is the following equation. 


u 0 - t n = Un - to = ^(mf - ml) (ml - m\) . 


(57) 
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Thus, for the reaction NN — > anything, we have 


(58) 

(59) 


uq t-ji , 

Utt to • 

We now consider the two cases of elastic and inelastic scattering. 
Elastic scattering (e.g. NN — » NN) 


For the case of m\ = m 2 = m 3 = m 4 = mjy, the threshold value of s, denoted St is 

s t = 4 m 2 N . (60) 

The equal masses case gives 

A12 = A 34 = s(s - s t ) = s(s - 4 m 2 N ) . (61) 

By equation (56), 

to = 0 (62) 

and 

t 7 r = st — s = 4m^- — s . (63) 

Clearly, to has a fixed value, but the value of t n depends on s. The minimum value of s is 
st = 4 m 2 N and so the minimum and maximum values of t n are 

C n (s->oo) = -00 (64) 

and 

C ax (s = s t ) = 0, (65) 

and so the range of t n is 

t-jr : —00 — > 0 . (66) 

Clearly, the singularity value t = m 2 is never encountered when t ranges between to and t n . 
Now, consider the singularity at u = m 2 , where s + t + u = m\ + + rn 2 + m 2 = 4 m? N = st- 

Thus, u = 4ra| — s — t = st — s — t. The extremum values of u can be obtained with equation 
(59). Thus, the extremum values of u are the (opposite) extremum values of t. Given that the 
singularity value t = m 2 is never encountered when t ranges between to and t w , it follows that 
the singularity value u = m 2 is never encountered when u ranges between uq and u n . 
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Inelastic scattering (e.g. NN — > IVA) 


For the case of m\ = m 2 = m 3 = 777,7V and 7714 = 777 a, the threshold value of s, denoted St is 

st = (m n + mA ) 2 (67) 

and 

A12 = s(s — ^m 2 N ) (same as elastic case) (68) 

and 


A 34 = (s - 777, 2 n - m 2 A ) 2 - 4m 2 N m 2 A . 

One uses equation (56) to obtain 

s{m\ + 3m 2 N - s) ± \J s{s - Am 2 N )\j{s - m 2 N - m 2 A ) 2 

The symbol ± means that a + sign is used for to and a — sign is used for t n . 
maximum values of to are 


to(U) - 2g 


(69) 


— Am 2 N m\ (70) 
The minimum and 


C n (s = s «) = rn N {m N - 777, a) 


(71) 


and 


giving the range of to as 


/.max 

r 0 


(s 


00) 


0 , 


to : 7Ti]\r{rnN — m/±) — > 0 . 


(72) 

(73) 


Note that if we had m n in place of m a, these results reduce to equation (62). The minimum 
and maximum values of are 

C in (s^oo) = -00 (74) 

and 

C ax (s = «t) = m N (m N - 777a) = C n / ( 75 ) 

giving the range of t n as 

t 7 r : —00 —7 7777 V (t 77- TV — 777 a) • ( 76) 


If we had 7777V in place of 7?7 a, these results reduce to equations (64), (65), and (66). Clearly, 
the singularity value t = t? 7 2 is never encountered when t ranges between to and t n . 

Now consider the singularity at 77 = m 2 . The extremum values of u can be obtained with 
equation (59). Thus, the extremum values of u are the (opposite) extremum values of t. Given 
that the singularity value t = m 2 is never encountered when t ranges between to and t n , then 
the singularity value u = m 2 is never encountered when u ranges between uq and u n . 
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2.4 Scattering via zero-mass exchange particle 

In this section, we consider scattering by means of an exchange particle with mass m = 0. By 
setting the coupling constant g = g^NN = a — 1/137, the exchange particle becomes the photon, 
except that in keeping with the scalar model, the photon is treated not as a spin 1 particle, 
but as a scalar. In the next section, both the elastic NN — > NN (Rutherford) and inelastic 
NN — > N A cases are examined numerically. The formalism developed in the preceding sections 
gives the well known differential cross section for Rutherford scattering, which is characterized 
by 9iZx = 92 Ax = 9, the exchange particle mass m x = 0, and in the general case m\ = m 3 
and m 2 = m .4 (i.e. elastic scattering). In our special case, for simplicity we make the further 
restriction mi = m 2 = Af, so that, in the cm frame, |pi| = | P 2 1 = |p 3 1 = | P 4 1 = |p cm | = |p|- 
Also, E\ = E 2 = = £4 = E. Thus, 

t=(p 4 -p 2 ) 2 = m\ + m\ — 2 E 4 E 2 + 2p 4 • p 2 

= 2 M 2 — 2 E 2 + 2p 2 cos 6 
= 2(— p 2 + p 2 cos 6) = — 2p 2 (l — cos 9) 

= ~Ys^ l ~ cos6 ^' ^ 

where, by equation (12), 

A = A (s, A/ 2 , AT 2 ) = s(s - 4M 2 ) = 4.sp 2 . (78) 

Also, 

u = (p:i - P 2) 2 = m| + — 2£ , 3 £ , 2 + 2p 3 • p 2 

= m\ + m 2 - 2E 3 E 2 - 2p 4 • p 2 

= -^(l + cos0). (79) 

Thus, by equation (9), the invariant amplitude, with m x = 0, is 


_ _ 2 ri 11 _ 2 r 1 1 

9 t + u 9 (P4 - P2) 2 + {P3 -P2) 2 . 

2.sg 2 1 1 

A 1 — cos 6 1 + cos 6 

4 sg 2 

A sin 2 9 
4 g 2 

(s — AM' 2 ) sin 2 9 

Equation (44) becomes 

= — — — |A4| 2 = — — , 

dQ cm 647t 2 s 647t 2 s ( s — AM 2 ) 2 sin 4 9 


( 80 ) 


(81) 


13 



(82) 


giving the Rutherford formula (with S = 1/2!) as 


da 

dQ 

err 


s 

647t 2 s 


\m\ 2 = 


87 t 2 s(s — 4M 2 ) 2 sin 4 9 


This agrees with Griffiths [7] (equation 6.54), as can be seen using s = 4 E 2 and p 2 = ^ = 

— 4^ - = j(s — 4M 2 ). The classical and quantum mechanical results for da /dVl exhibit several 
differences (compare Griffiths [7] equation 6.54 with p. 194). The classical result for Rutherford 
scattering typically encountered (e.g. Griffiths [7] (p. 194), or Goldstein [11] (equation 3.102) 
depends on 9/2 rather than 9. This difference of a factor of 2 is due to the fact that the classical 
calculation is usually performed in the lab (target) frame, whereas the present calculation is 
performed in the cm frame. In the non-relativistic limit, for the case of all particles having the 
same mass, we have precisely 9i a b = 9 cm / 2. See Marion [12] (equation 9.71). 

The classical and quantum mechanical results for da/dVL also differ by an overall factor of 
2. This well known difference arises from the quantum mechanical description of a particle as a 
wave packet, versus the classical description of particles as having finite extent (see Bohrn [13] 
sections 17 and 49 for a discussion). As expected, the Rutherford formula diverges at s = 0 
and s = 4M 2 , just as the formula of Griffiths [7] which diverges at E = 0 or p 2 = 0. Since 
s = (p\ +P2) 2 ) the value of s in the cm frame is s = (E\ + E 2 ) 2 , which has a minimum value at 
s = 4 M 2 for mi = m 2 = M. Thus, the singularity at s = 4 M 2 is the same as the singularity 
at p 2 = 0. As for the case s = 0, this value is never reached because the minimum value of s 
is 4 M 2 . Similarly, the value E = 0 is never reached because the minimum value of E is M. As 
mentioned in the previous section, the equation for the angular distribution is much simpler if 
Mandelstam variables are used instead of angle, 


da Sg 4 1 1 2 

dQ cm 647 r 2 s t u 

with u = 4 M 2 — s — t. 


(83) 


2.5 Total cross section 


In this subsection the previous differential cross sections will be integrated to form total cross 
sections. The total cross sections are obtained by integrating [10] 


or 



a = 



da 

du 


(84) 


(85) 


The expressions for the direct, exchange, and interference t distributions were given in equations 
(16), (17) and (18). The total cross sections can be obtained with straightforward analytic 
integration. The direct cross section is 





K 

(t — m 2 ) 2 


= -K 


1 


t — m 2 


( 86 ) 
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giving 


a d =K 


f, r — m 2 to — m 2 


(87) 


The exchange cross section is 


fio 


= 


dt 


K 


u («-m 2 ) 2 


fio 


dt 


I< 


rto 


dt 


I< 


l tn (Tim 2 — s — t — m 2 ) 2 J tn (i + s-Sm 2 + m 2 ) 2 

1 


—K 

K 

K 


" to 

- t-TT 


1 


t + s — Tm 2 + m 2 

1 

+ s — Tm 2 + m 2 to + s — Tm 2 + m 2 

1 1 


t n r — m 2 to — m 2 


giving 


— Gd, ■ 


( 88 ) 


(89) 


Thus, the direct and exchange total cross sections are equal. Note that this result is more easily 
obtained by direct integration of the u variable. Since du = —dt, and using (59), we have 


oy = 


[ u * J K 

L (« - "» 2 ) 2 

fto dt K 


J u ( u - m 2 ) 2 ■ 

= CTd 

Now consider integration of the interference term given in equation (18). 

d(Ji 2 K 2 1\ 

dt ( t — m 2 )(u — m 2 ) ( t — m 2 )(—t + Tm 2 — s — m 2 ) 

2 K 

—t 2 + t(Tm 2 — s) + m 2 (s — Tm 2 + m 2 ) 

The needed integral is of the form [14] 


(90) 


(91) 


dx _ 2 

ax 2 + bx + c a/4 ac - b 2 


tan 


-l 


2 ax + b 
V^ac — b 2 


1 ( 2 ax + b — \Jb 2 — 4ac 


Tog 


\Jb 2 — 4ac, y2 ax + b + \Jb 2 — 4ac^ 

-1 


ax + ^b 


for b 2 — 4 ac < 0, 
for b 2 — 4 ac > 0, 
for b 2 — 4 ac = 0. (92) 
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Note that if b = 0, the result of equation (92) reduces to the correct form for both cases of the 
tan -1 or log. However, if a = 0 or c = 0, the result of equation (92) does not reduce to the 
correct form. In fact, for a = 0, we have [14] 


and for c = 0 we have [14] 

For our case, we have 
a = — 1 


dx 


bx + c 


dx 


- log (bx + c ) 
b 


1 , I x 

r lo § 


ax 2 + bx b \ax + b 


b = Em 2 — s 


c = m 2 (s — Em 2 + m 2 ) , 


(93) 


(94) 


(95) 


and we never come across the situation where a = 0 or c = 0. The minimum value of c occurs 
at s = Sthreshoid = ( m 3 + m 4 ) 2 . Thus, 

Cmin = m 2 [(m 3 + m 4 ) 2 — m 2 — m 2 — m 2 — m 2 + m 2 ] 

= m 2 — (mi — m 2) 2 + 2(m-3m 4 — m-im^) , (96) 

but since we always have m 3 ?n 4 > m±m 2 , then 

Cmin > rn 2 - (mi - m 2 ) 2 , (97) 


giving 


m = m n > |mi — m 2 | => c > 0. (98) 

Thus, we can conceivably run across the case where c = 0. Now we can have b = 0 at s = Em 2 
but that is no problem because, in that case, the integral above reduces to the correct form. 

Given that the result depends on the value of b 2 — 4 ac, we need to evaluate this quantity for the 

expression in equation (91), with the result that 

b 2 — 4 ac = (Em 2 — s) 2 + 4 m 2 (s — Em 2 + m 2 ) . (99) 

Define s = s — Em 2 + m 2 so that s — Em 2 = s — ?n 2 , to give 

b 2 — 4ac = (s — m 2 ) 2 + 4m 2 s = s 2 + 2m 2 s + m 4 , 

= (s + m 2 ) 2 (100) 

which finally gives 

b 2 — 4 ac = (s — Em 2 + 2 m 2 ) 2 , (101) 

which, being a square, always gives 

b 2 - 4ac > 0 . (102) 


16 



Following a similar analysis for c above we have 


(i b 2 — 4 ac) min = [(m3 + m ^) 2 — m 2 — m\ — m 2 — m 2 + 2m 2 ] 2 

= [2m 2 — (mi — m2) 2 + 2(777,3771.4 — mim2)] 2 , 

but since we always have m3 7774 > mi m2, then 


giving 


(b 2 - 4 ac) min > [ 2 m 2 - ( 777-1 - m 2 ) 2 ] 2 , 


\plm> I mi — 7772 ! => fr 2 — 4 ac > 0. 


Thus, we can conceivably run across the case where b 2 — 4ac = 0. The following analysis assumes 
that c / 0 and b 2 — 4 ac 7 ^ 0. The total interference cross section is given by 


a j = dt 


t —t 2 + i(Sm 2 — s) + m 2 (s — Em 2 + m 2 ) 


2 K 

s — Em 2 + 2 m 2 


t + s — Em 2 + m 2 to 
t — m 2 , 


which becomes 


2K (t n — m 2 )(to + s — E 777 2 + m 2 ) 

s — Em 2 + 2m 2 (to — m 2 )(t,r + s — E 777. 2 + m 2 ) 


Elastic scattering (e.g. NN — > A r 7V) 


For the reaction iViV —* NN, we obtain 


to — 0 , 

U = 4 m 2 N - s , 


which, upon substitution into equation (87), gives 


0(1 ^ ?77 2 S — 4777-^r + 777- 2 

VttTVJV ( 1 


32-7ts(s — 4mi) \?77, 2 s — 4mlb + m 2 


which simplifies to 


ov = 


9ttNN 


327 t m 2 s(s — ^m 2 N + m 2 ) 
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Note that the minimum value of s is st = ^thr esh old = (m 3 + 7714) 2 = 4 m? N and so the cross section 
has no singularities. The largest value of a occurs at s = St- 


&d(s -> s t ) 


9nNN _ 9 n NN 
327T m^st 1287T m^vn? N 

1 / 9-kNN 

&Trm 2 N \ 2 iti 7V 

S ( g-n-NN 

4:Trm 2 N \ 2m n 



The other important limit is 


cr d (s -*■ 00 ) = 0. 


( 112 ) 


(113) 


Recall (89), which showed that the direct and exchange total cross sections are the same. The 
sum of direct and exchange maximum values is 


^max 1 ^max 
°d + ® e 


2 S / 9ttNN 

4itm 2 N \ 2 m n 

1 ( 9nNN 

47T m 2 N \ 2m n 


4 


4 


The interference cross section is given by substitution into equation (107), 


(114) 


which is 




2 K 

s — 4m^ + 2 m 2 
4 K 

s — 4 m? N + 2 m 2 


log 

log 


(4m^ — s — m 2 )(s — 4 m 2 N + m 2 ) 
_—m 2 {4m 2 N — s + s — 4 m 2 N + m 2 ) 
( s- 4 m 2 N + m 2 \ 

V rnl ) 


0 * = 


9nNN 


8irs(s — 4 m 2 N )(s — 4 m 2 N + 2m 2 ) 


log 


s - 4 m 2 N + 


m 


m* 


= \a d m 


s - 4 m 2 N + 


m 


1 


s — 4 m 2 N + 2 m 2 s — 4 m\ 


log 


s - 4 m 2 N + 


m 


N 


m 


(115) 


(116) 


This cross section has a potential threshold singularity at s = st = 4 m 2 N , which comes from 

4 

the term K = K(s) = j£x^9i 3x 92ix = 32 ns{s-l mf,) The sin g ularit y was canceled out in our 
previous expression for a d in equations (110) and (111). Clearly, we need to evaluate the term 
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log 3 ^ 2 ^* near s = st- Consider the following Taylor series expansion [14] (p. Ill) 


(« - st) 


log 


s — st + m z 


m * 


1 


(» - s t) 

2 log(l T x) 
m z x 

1 


log(^ + l 


to 


with x = 


s- s t 
m 2 


m 2 x 


2 S 4 
X X° X 


x ~y + y _ T + " 


If x x 2 x 6 

n? I ~ 2 + T“T 


which gives 


1 


iinr 

s^st (s - s t ) 


log 


S — St + 77T 


TO. 


1 


TO 


2 ' 


(117) 


(118) 


Substituting into the above expression gives the interference cross section at threshold s = st, 

= (Ti(s — > St) = 2<J U d 


TO; 


r max 

T d 


^max | ^max 

= CTrf + <TO 


(119) 


which is an interesting result because it shows that in the low energy limit, the interference term 
dominates. The other important limit is 


&i{s — > oo) = 0 . 

Inelastic scattering (e.g. NN — > IV A) 

For the NN N A reaction, we obtain from equation (56) 


to(U) - 2s 


s{m\ + 3 m 2 N — s) ± \J s(s — 4to^) \J (m\ + m 2 N — s) 2 — 4 m 2 N m. 


(120) 


( 121 ) 


Substituting gives 


crd = 


91nn91na \J(s- m\ - to 2 ,) 2 - 4 m 2 A m 2 N 


and 


TO; = 


167T \J s(s — 4 m 2 N ) s 2 m 2 + sm 2 (m 2 — m 2 A — 3 m 2 N ) + ( m 3 N — m^mjv) 2 


2 2 
9-kNnS-kNA 


47ts(s — 4m^)(s + 2 to 2 — m 2 A — 2>m? N ) 
x log 


s(s + 2 to 2 — to^ — 2>m 2 N ) + J s(s — 4 m? N ) J ( s — m 2 A — m 2 N ) 2 — 4 m A m 


(122) 


N 


s(s + 2 to . 2 — m A — 3 m 2 N ) — yj s(s — 4 m 2 N ) yj (s — m 2 A — m 2 N ) 2 — 4 m 2 A m 2 N 


(123) 
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Note that if one replaces the A mass with the nucleon mass, these results reduce to the elastic 
formulas obtained previously (except for the statistical factor S). Also note the following limits. 


* ^threshold) 

* ^threshold) 


cr d (s — > oo) = 0 , 
CJi(s — > oo) = 0 . 


(124) 

(125) 


2.5.1 Asymptotic region 

We now consider a theoretical check of our cross section result. The Froissart theorem [15] (p. 
344) states that in the high energy limit s — > oo, the total cross section of an arbitrary collision 
process is bounded by 

cr tot < Constant x [logs] 2 . (126) 

For the elastic pp —> pp and inelastic processes pp — » pA + , allowing s to become large leads to 
equations (113), (120), (124) and (125), which shows that the total cross sections predicted by 
the scalar OPE model go to zero in the high energy limit, in agreement with the theorem. 

3 Calculations for 2 - body final states 

The previous section developed the basic theory for calculating 2 - body final state cross sections. 
In the present section, these cross section formulas are illustrated with numerical results and 
plots for the scattering processes pp — ► pp and pp —> pA, and the decay A — ► pn°. 

3.1 Coupling constants 

In this subsection, some comments will be made about coupling constants. In quantum field 
theories, the coupling “constants” are in fact not constants, but depend on the momentum p of 
the interaction, thus (j-kNN = g-nNN^p) [7] (p- 62), [16] (pp. 146 - 151). From the perspective 
of the quark model, g^NN hides the underlying quark-gluon interactions, whose strengths are 
a function of the momentum of the quark-gluon interactions. The low energy pion-nucleon 
coupling constant g W NN has been experimentally determined to be [17] 

&nn/ 4?t « 14 , (127) 

for the coupling of neutral pions to protons and neutrons. The coupling of charged pions is a 
factor \/2 stronger [18] (p. 184); [19] (p. 219). In quantum field theories where spin is taken 
into account, the coupling constants are dimensionless. However, in the scalar OPE theory the 
coupling constants must have dimensions of GeV [7] (p. 202). 

3.2 Calculations for scattering involving a massive exchange particle 

The exchange particle is always either massive or massless, and quite different behaviors in cross 
sections arise as a result. In this subsection, we give results for massive particle exchange. 


20 



3.2.1 Minimum and maximum values of t 


Consider the reactions pp — » pp and pp — > pA + . We ignore the particle spins and assume that 
the exchange particle is a neutral pion in both cases. This is consistent with the use of a scalar 
theory, in that other pion isospin states are not included. Also, the use of a single pion exchange, 
as opposed to multiple exchanges, is an approximation expected to be valid at low energy. At 
higher energy, multiple exchanges can be introduced, or parameters, such as coupling constants, 
can be adjusted to fit data. Using the following particle masses m p = 938.27 MeV, m x = m^o = 
134.98 MeV, m&+ = 1232 MeV, we wish to calculate the value of in units of nrb/GeV 2 for 
both reactions at s = 2s t hreshoid (twice the threshold to produce the A) for the maximum and 
minimum values of t given by (56) 


to(^7r) 


1 

4 s 


( m\ — m% 



(128) 


The notation t n is bound to be a little confusing given the the particle we are talking about is 
the pion with symbol ir. The symbol refers to the value of the t variable at the angle of 7r radi- 
ans. We obtain the following results. For the reaction pp — > pp, we calculate the following results. 


da 

dt 

da 

dt 


■(» 

(« 


•^threshold 
U (^threshold) 
Ur ( 2 St hreshold ) 


^SthresholdT — U)) 


^Sthreshold) f — ^7r) 


3.52 GeV 2 , 

0 GeV 2 , 
-3.52 GeV 2 , 


14700 

14700 


mb 

GeV 2 

mb 

GeV 2 


For the reaction pp — > pA + , we calculate the following results. 


(129) 

(130) 

(131) 

(132) 

(133) 


da 

dt 

da 

dt 


■(» 

(« 


•^threshold 


to ( 2 ^threshold ) 

Ur (^threshold) 


^Sthresholdi t — t o) 


^Sthreshold) t — U ) 


4.71 GeV 2 , 
-0.00723 GeV 2 , 
-5.25 GeV 2 , 


6720 

6720 


mb 

GeV 

mb 

GeV 


2 > 


2 • 


(134) 

(135) 

(136) 

(137) 

(138) 


The cross sections are the same at the values of to and Ur because the center of momentum 
frame does not distinguish forward and backward scattering. 


3.2.2 Singularities 

The invariant amplitude M. in equation (9) has singularities at t = m 2 and u = m 2 . These 
singularities are not encountered for the reactions considered above. To show this, we evaluate 
t at the singularities and compare to the values of to and t n for the case of s = 2s t h res hoid- 
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Results for the reaction pp — > pp are 



^threshold — 3.5214 GeV 

(139) 


^o(2sthreshold) = 0 GeV^ 

(140) 


G ( 2 ^threshold ) = —3.5214 GeV 2 . 

(141) 

The t singularity occurs at 

t = ml = 0.0182196 GeV 2 . 

(142) 

The u singularity occurs at 



7t = 777 2 + 77 t| + 777§ + 777-4 — m x ~ S ~ t = —3.53962 GeV 2 . 

(143) 


For these maximum and minimum values of t (to and G)> the singularities are never encountered. 
Results for the reaction pp — > pA + are 

^threshold = 4.71007 GeV 2 , 
to(2sthreshoid) = —0.00722818 GeV 2 , 

G(2sthreshoid) = —5.25404 GeV 2 . 

The t singularity occurs at 

t = m 2 = 0.0182196 GeV 2 . (147) 

The u singularity occurs at 

u = m\ + m\ + m\ + m 2 — m 2 — s — t = —5.27949 GeV 2 . (148) 

For these maximum and minimum values of t (i.e. to and t n ), the singularities are never en- 
countered. 

3.3 Invariant distribution da/dt 

The Lorentz invariant differential cross section da/dt is very useful because it represents the 
complicated, non-invariant angular distribution in a very simple way. This subsection develops 
explicit formulas for this invariant cross section. Equation (11) is used to calculate the invariant 
distribution da/dt. The total cross section a is determined by integrating da/dt over the range 
t = [io,G]> with to and t n given by (56). The exchange particle is the 7 r°. By setting the 
mass 777.4 = 7/7 a , the mass of the A particle, we generate results for the reaction pp — > pA + . 
For the reaction pp — > pp, results are obtained by setting 777.4 = t n p , the mass of the proton. 
Figure 4 plots da/dt at s = 2s t h re shoid- The shape of the plot agrees with Reference [20] (figure 
4.8). The t and u singularities are prominent at the left and right extremes of the plot. The 
direct and exchange terms are plotted separately in figures 5 and 6. The total cross section at 
2sthreshoid is a(2st) = 40.2605 mb. The total cross section a is plotted in figure 7 as a function 
of s, and has a maximum of ~ 90 mb. Since the range of integration lies between the t and 


(144) 

(145) 

(146) 
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u singularities, the Mathematica function NIntegrate is able to integrate the curve as a single 
piece, unlike the case involving da/dClsiab (see section 3.6), where the singularity at 0 pea k requires 
piecewise integration, a goes to zero at Sthreshoid- As s increases, a rises to a peak and then 
slowly diminishes with increasing s. These characteristics of the shape of the curve predicted 
by the scalar ABC theory agree with experiment as seen in reference [21] (figure 40.11). 


p+p->p + A + (S = 2 Sthreshoid) 



Figure 4: Differential cross section da/dt. 


Direct term: p+p->p+A + (s = 2 Sthreshoid) 



Figure 5: Differential cross section da/dt direct term. 
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Exchange term: p+p->p+A + (s = 2 s t hreshoid) 



Figure 6: Differential cross section da/dt exchange term. 

p+p->p+A + 



s (GeV 2 ) 


Figure 7: pp —> pA + total cross section. 


3.4 Angular distribution d<j/dfl c and o in cm frame 

The angular distribution in the cm frame will be developed in this subsection. For the 2 - body 
reaction 1 + 2 — » 3 + 4, the angular distribution in the cm frame is given by equation (44), 


da 




cm 



— \J ^12^34 

47TS 


da 
dt ’ 


(149) 
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with 


t = ml + m\ + 

Ah 


- \J (A12 + 4sm|)(A 34 + 4sm|) + \J A12A34 cos 0 2 4c 


(150) 


See equation (45). In the cm frame, the angular distributions of particles 3 and 4 are equal [10], 

dcr/dQzc = da/dQ^c . (151) 


Since both forms of the angular distribution above are functions of t, either may be used in 
conjunction with the expression f(024c) to give da/dtl^c as a function of 024c- The exchange 
particle is the 7r°. By setting the mass 777.4 = 777, a , the mass of the A particle, we generate results 
for the reaction pp — > pA + . For the reaction pp — > pp, results are obtained by setting 7714 = m p , 
the mass of the proton. In figure 8, the plot of da /d£l\ c at s = 2sthreshoid exhibits a similar shape 
and the same symmetry as the plot of da/dt, noteworthy since one is a function of 6 while the 
other is a function of t. Evaluating da/dQ^c at 6 = 0 and 6 = n gives 


da 

dVL^c 

da 

d&4C 


:(» 

(» 


2s t ,0 = O) 
2 s t , 6 = 7 r) 


mb 

318.376 — , 
sr 

mb 

318.376 — . 
sr 


(152) 

(153) 


The total cross section a is found by numerically integrating da/dVL^c over all possible final 
state angles 0 = [0, tt] . At s = 2 st, a(s = 2 s*) = 40.2605 mb, which matches the previous 
result found by integrating da/dt. Figure 9 shows the total cross section cr(s) over the range 
s = [^threshold) 2s t h res hoid] • The plot is identical to the plot resulting from the integration of 
da/dt , as expected since total cross section is invariant. 


p + p >p + A ( S — 2 Sth res hold) 



Figure 8: Differential cross section da/dVt 4 C . 
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p+p->p+A + 



s (GeV 2 ) 


Figure 9: Total cross section a. 


3.5 Spectral distribution da/dE 3 \ in the lab frame 

In this subsection, the spectral distribution da / dE 3 i is derived from the invariant distribution 
da/dt by expressing the invariants t and u in terms of the energy E and kinetic energy T, 
respectively (see Section 2.3.3). The formula is [10] 


da da da 

dEi = dEu = + 1712 ~di 


(154) 


with t given by equation (37). The exchange particle is the i r°. By setting the mass 777,4 = m a , 
the mass of the A particle, we generate results for the reaction pp — * pA + . The results for the 
reaction pp — > pp may be obtained by setting 777,4 = m p , the mass of the proton. The differential 
cross section da/dE 4 ; is plotted at s = 2 ,s threshold • See figure 10. Note that in reference [10], it 
was shown that dE 3 i = —dE^i and dT 3 i = — dT, 4 ;. Therefore, plotting da/dE 4 ; is equivalent to 
plotting da/dE 3 i. The plot range is found by inverting the formulas for t and u as functions 
of energy and evaluating them at to and t n . The plots have the characteristic horseshoe shape 
of the invariant distribution (which is plotted as a function of t) . The total cross section is 
determined by numerically integrating the spectral distribution, and is plotted as a function of 
s over the range ^threshold to 2s t hreshoid- At s = 2sthreshoid, the total cross section 


a(2s t ) = 40.2605 mb 


(155) 


is identical to the value determined using the invariant distribution. Note that the two-body 
lab frame spectral distribution in equation (154) does not possess any singularity. Also, recall 
that a two-body cm frame spectral distribution is impossible to form because the cm energy of 
the final state particles are fixed [10]. In figures 4 and 5 of reference [22], Murphy, Dernier and 
Rarnaty present cm and lab frame spectral distributions for the reaction p + a — * 7r° + X. These 
are very useful data for future comparisons to theoretical models. Both spectral distributions 
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are finite and well behaved over the entire range of pion energy. However, the results above do 
not apply to this data because the reaction p + a — > ir° + X is at least a 3-body final state. 


p + p->p+A + (S = 2 Sthreshold) 



E 4 lab [GeV] 


Figure 10: Differential cross section da/dE 4 ;. 


3.6 Angular distribution da/dVt 3 i in lab frame 

The angular distribution is required in the lab frame for space radiation studies. This is devel- 
oped in the present subsection. The angular distribution da/dQ 3 i nray be expressed in terms 
of the cm distribution da/dVL 3c . This is achieved by expressing t in terms of £ 3 ; , and £ 3 ; in 
terms of 9 3 ;. The exchange particle is the 7 r°. By setting the mass 1714 = uia, the mass of the A 
particle, we generate results for the reaction pp — ► pA + . The results for the reaction pp — > pp 
may be obtained by setting 1714 = rn p , the mass of the proton. For the reaction pp — > pA + with 
particle 3 the A + particle, the energy E 3 i(9 3 i) is a double valued function, except at # pea k where 
the two energy roots are equal. In figure 11, the two £ 3 / roots are plotted piecewise together and 
join smoothly at $ 3 / pea k forming a single curve. E 3 i root i(9) forms the upper part of the curve, 
while E 3 i root2 (9) forms the lower part of the curve. Since £ 3 h-ooti(^pe a k) = £ 3 / 100 * 2 (# pea k), ^ pea k 
is found by setting the square root term in the numerator of (51) equal to zero and solving for 
9. The slope of £ 3 ; goes to infinity at 0 pea k- This causes a singularity in da/dQ 3 i at 0 pea k- The 
presence of the singularity is demonstrated by expressing da/dQ 3 i in terms of dE/d9, 

da da da 1 dE 

d£l 3 i d{cos9)dcj) dE (— sin 9) drf> d9 

da/dE 3 i is continuous over the allowed range of E- 3 i (and corresponding range of 9 ), and 
l/sm# pea k is finite at 0 pea k, whereas dE/d9 — * 00 at ^peak- Physically, a narrow detector 
situated at an angle 9 < 0 pea k will detect particles at two energies. The double valued nature 
of E 3 i(9) necessitates that da/dQ 3 i be evaluated in two parts, one part as a function of £ 3 z r ooti> 
the other of E 3 i TOOt2 . da/dQ 3 i is the sum of the two parts. This process represents the adding 
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together of the contributions to the differential cross section of distinct states, which are physi- 
cally distinguished by their different energies. The cross sections are added together only after 
their respective amplitudes are squared. In figure 12, the plot of da/dQ^l shows the singularity 
at 03 /peak — -744. The partial cross section Aa is defined to be the cross section integrated over 
a finite small angular range. In figure 13, the detector plot shows the partial cross section Act 
seen by a narrow detector of angular width Ad = .01745, positioned at angles 631 ranging from 0 
to $ 3 i P eak' The plot reveals an increase in the partial cross section near 03 i pea k- The relativistic 
kinematics predict a bump in the lab angular distribution in the vicinity of 03 z pea k- The plots 
of t and u as functions of 631 support the discussion about direct and exchange terms in section 
2.3.1. For a given angle 631 and energy E 31 , there corresponds two unique values of t, and of u. 


E31 vs 03 



03 (rad) 

Figure 11: Energy £ 31 ( 631 ). Upper and lower parts of the curve correspond to E 31 roots 1 and 
2, respectively. 


3.6.1 Relation between cm and lab angles 

In the cm frame, particles 3 and 4 exit from the collision in opposite directions, at angles ranging 
from 0 to 27 t. Angles are measured with respect to the collision axis (particle 1 moves parallel to 
the collision axis in both frames, particle 2 is at rest in the lab frame but moves parallel to the 
collision axis in the cm frame). Due to symmetry about the collision axis, only angles between 
0 and 7T need be considered. When viewed from the lab frame, the exit angles of particle 3 
and 4 differ from the cm angles as shown in figure 14, where cm angles are transformed into 
corresponding lab angles. The difference in angles between the cm and lab frames is due to the 
relative velocity of the two frames, which affects the component of particle momentum along the 
collision axis, but not the component transverse to the axis. Also, see figures 15 and 16 which 
show the Mandelstam variables as a function of lab angle. 

For the reaction p + p — ► p + A + , with particles labeled 1 + 2 — > 3 + 4, the velocities of 
initial state particles 1 and 2 in the cm frame are equal and opposite, and the velocity of particle 


28 




p + p->p+A + (S = 2 Sthreshold) 



Figure 12: Differential cross section da/dQ^i. 


5 

4 

dela(mb) ^ 
2 

1 

0 


p+p->p + A (S = 2 Sthreshold) 



Figure 13: Detector partial cross section. 


2 is equal to the velocity in the lab frame. In the cm frame, the momenta of the final state 
particles 3 and 4 are equal and opposite, but the velocities differ, thus the speed of the lab frame 
may be less than or greater than the speed of a final state particle. For the case of particles 3 
and 4 exiting along the collision axis (with angles of 0 and 7 r), this means that in transforming 
angles from cm to lab frame, the direction of a particle may be reversed (the particle with angle 
9 cm = 0 in the cm frame may have angle 8i a b = 7r in the lab frame). For all other angles, the 
transform from the cm frame to the lab frame reduces the angles, that is 9i a b < 0 cm . The effect 
of transforming the angles on the angular distributions is seen by comparing figure 8 with 12. 
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1 031 1 and |041j versus 0c 



Figure 14: Angles 0 3 i ab and Q^ ab versus 9 cm . 


t as function of 03 



Figure 15: Mandelstam variable t as a function of 6^. 

3.7 Calculations for zero-mass exchange particle (Rutherford scattering) 

Previous sections were concerned with massive particle exchnage. In this subsection, it is shown 
how to handle massless particle exchange. In Rutherford scattering, for the case of initial and 
final state particles having the same mass, the exchange particle has a mass m = 0 , and the 
coupling is electromagnetic with coupling constant (J^nn = a (in place of the pion coupling 
constant g-n-NN)- The exchange particle is the photon 7 with mass m = 0. For pp —> pA 
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u as function of 03 



Figure 16: Mandelstam variable u as a function of O 31 . 


(non- Rutherford), evaluating da/dt at to and G at s = 2 s t hreshoid gives 

^rr{s = 2s t hreshold,i = t 0 ) = 7.59 P \ , (157) 

at Gev 

7 , (s 2sthreshold) 1" G) 7.59 — _ o • (158) 

at GeV 

For pp — > pp (Rutherford), these cross sections are infinite due to the t and u singularities 
at t = to and t = t n . The physical explanation of this is that scattering always takes place 
when there is a potential of infinite range. Figures 17 and 18 plot the differential cross sections 
da/dt and da/dQ cm for pp — > pA + . Figures 19 and 20 plot the differential cross sections for 
pp — > pp. Figure 21 plots the total cross section for pp —> pA + (non- Rutherford) . For pp — > pp 
(Rutherford), the total cross section is infinite, in agreement with the classical result. 


3.8 Comparison of cross sections 

In this subsection, different forms of the cross section will be compared. The relative strengths 
of interaction processes pp — > pp and pp — >■ pA + , both for massive and zero-mass exchange 
particles, are shown by comparing plots of the cross sections. For Rutherford scattering, the 
total cross section is infinite and cannot be plotted. Figure 22 compares plots of total cross 
sections for the elastic pp — > pp and inelastic pp — > pA + processes involving a massive (pion) 
exchange particle, over the energy range Sthreshoid to 2s t h re shoid- The elastic process dominates 
at low energies. Figures 23 and 24 plot the invariant cross sections of the elastic process pp — ► pp 
for massive (pion) and zero-mass (photon) exchange particle. Judging by the vertical axis scales, 
the strength of the pion process is many orders of magnitude greater than the photon process 
over the valid range of variable t. The plots are limited in height, and do not reveal the full story. 
The Rutherford invariant cross section goes to infinity at the limits to an d t n , whereas the pion 
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p + p->p+A + ( S = 2 S t hreshold) 



Figure 17: Invariant differential cross section (non-Rutherford). 

p + p->p + A + (S = 2 Sthreshold) 



Figure 18: da/d£l cm Differential cross section (non- Rutherford) . 


invariant cross section is finite at the limits, da/dt = 37.7 rnb/GeV 2 . Similar statements apply 
to the cm frame angular distributions, which are closely related to the t distributions since t 
varies linearly with cos 6. See equation (45). For most of the allowed range of angles 0 to 7 r, the 
angular distribution of the pion process is many orders of magnitude greater than the photon 
process. However, at angles 0 and 7r, the Rutherford angular distribution goes to infinity. This 
is apparent by equation (81). 
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p + p->p + p (S = 2 Sthreshold) 



Figure 19: Invariant differential cross section (Rutherford). 

p+p->p+p (S = 2 Sthreshold) 



Figure 20: da/dQ cm Differential cross section (Rutherford) 


4 Scattering for 3 - body final states 

In this section, the formalism for describing 3 - body final state reactions will be presented. 
Cross sections will be calculated and special attention will be paid to the intermediate resonance 
formation and decay. Approximate techniques for handling resonance decay will be introduced. 
These results are important because many space radiation reactions include three or more bodies 
in the final state. The reaction 


NN -> NNir , 


(159) 
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p+p->p+A" 



s (GeV 2 ) 

Figure 21: Total cross section (non-Rutherford). 

7T 7 T . 

pp -> pp vs pp -» pA 



s (GeV 2 ) 


Figure 22: Comparison of elastic and inelastic total cross sections. The steadily falling cross 
section is elastic, and the cross section which rises and then falls is inelastic. 


proceeds through A resonance formation 

NN — > 1VA — > NNn . (160) 

The Feynman diagram for this reaction is shown in figure 25. Thus, the reaction proceeds in 
two steps, namely 


NN -»• N A , 


(161) 
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7T 

P+P -> p+p 



t (GeV 2 ) 


Figure 23: Invariant distribution for pp — > pp by pion exchange. 


da 

dt 


( 


mb 

GeV 2 


p+p -» p+p 



t (GeV 2 ) 


Figure 24: Invariant distribution for pp — > pp by photon exchange. 


and 


A — > Ntt 


(162) 


4.1 Amplitude for A resonance formation and decay 

An essential piece to the calculation of 3 - body final state cross sections, is the amplitude for 
resonance formation and decay. This is considered in the present subsection. The amplitude for 
the reaction in equation (159) is written in terms of the amplitudes in equations (161) and (162). 
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Figure 25: The Feynman diagram for the reaction NN — > N A — > NNir. Note the labeling of 
the particle numbers. The initial state nucleons are labeled as particles 1 and 2. The final state 
nucleons are labeled as particles 3 and 5. The intermediate A state is labeled as particle 4, and 
the final produced ir is labeled as particle 6. 

Consider the complete reaction NN — > N A — > NNtt. The Feynman diagram is illustrated in 
figure 25. Using the Feynman rules for ABC (scalar) theory gives the amplitude 

iM (27r) 4 J 4 (pi + p-2 - p 3 - P5 ~ Pg) 

= {-ig-KNN)(-ignNA) 2 {^) 4 S 4 (p 1 -P3- q n ) (2ir) 4 8 4 (q^ + p 2 - 9A)(2vr) 4 5 4 (g A ~ Ph ~ Pg) 

1 f ,4 I f ,4 i i 

X (2vr) 4 J Q7T (2vr) 4 J C q 2 - m 2 q\ - m\ 

= ( ig-rrNN ) ( A ) 2 (27t) 4 (5 4 (pi ~ p 3 + P2 ~ <?a) (27t) 4 <5 4 (^A “ P5 ~ Pg) 

1 f ,4 * E 

X (2vr) 4 J QA {p\ - p 3 ) 2 - m 2 q\ - m\ 

= (-igTTNN){-igTvNA) 2 (2ir) 4 S 4 (pi -p 3 +P2 -P5 - Pg) 
i % 

X (p\ - PgY - ml (p 5 - p 6 ) 2 - m\ ’ ( ' 163 ^ 

which reduces to 


iM( 1 + 2— >3 + 5 + 6 ) = iM(NN — > NNn) 

= {—igTrNN){ — ig-KNA ) 2 -y y -y y , 

TA~ m A 


(with q n = pi — p 3 and q A = Pg ~ Pg)- 
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Now consider the simpler reaction NN — > N A. The Feynman diagram is illustrated in figure 
26. The amplitude is 


■iM{2ir) 4 5 4 (pi +P2-P3- Pa) 

= ( ip-n-NN ) ( ( 27 t ) 4 d 4 (p\ - P3 ~ Qn) (2Tr) 4 6 4 (q 7T + P 2 - Pa) 

X (2tt) 4 J d Qn ql- ml 

= ( ^9ttNN ) ( ^ttN a) (2TT ) 4 <J 4 (pi -P3+P2 -P 4 ) 7 To T , 

(pi ~P3r - 


which reduces to 


(165) 


iM( 1 + 2 — > 3 + 4) = iM(NN — > N A)dLips 

'i 

= ( ’i , 9nNN)( ^SttN a) -n T , (with q n = pi — p 3 ). (166) 

K - ml 



Figure 26: The Feynman diagram for the reaction NN — > N A. 

Now consider the decay A — > Ntt. The Feynman diagram is illustrated in figure 27. The 
amplitude is 

zAl(27r) 4 (i 4 (p4 - p 5 - p e ) = (-ig 7r NA)(2TT) 4 5 4 (p 4l - p 5 - p 6 ) , (167) 

which reduces to 

iM( A — > Ntt) = —ig n NA • (168) 

Combining the above equations (164), (166) and (168) gives 

iM(NN -> NNtt) = iM(A -> Ntt) ^ iM(NN -»• IV A) . (169) 

Qa ~ m A 
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Ps=Pn 



Figure 27: The Feynman diagram for the decay A — ► Nit. 


Thus, 


1-^^ 12 — >356 1 2 — | -A/f 4 — > 56 1 2 W~2 2\2 1-^ 12^34 1 2 • (170) 

V PZ ~ m 4 ) 

The equations (169) and (170) express the amplitude of the complete process in terms of the 
sub-processes. 

4.2 Cross section for resonance formation and decay 

The amplitude of the preceding section is now used to determine the cross section. In the 
following discussion, all statistical factors are 5 = 1. Note that for NN —> NNir, 5 = 1/2. 
From reference [10], the cross sections are 

da(l + 2 -► 3 + 5 + 6) = -^| Afi 2 ^356| 2 (2vr) 4 d4> 3 (pi + p 2 ]P3,P5,P6) , (171) 

and 

dcr(l + 2 -4 3 + 4) = T ^|AIi2^34| 2 (27r) 4 d<I>2(pi +P 2 ;P 3 ,Pp • (172) 

4 r 

The phase space factor cf<h n is given in reference [10]. From reference [10], the decay width is 
cT(4 — > 5 + 6) = — \M^56\ 2 (2ir) 4 d$2(P4;P5,P6) • (173) 

27B4 

From reference [10], the phase space recurrence relation is 

d$ 3 (P;p l ,p 2 ,P3) = d$2(P-,q,P3)d®2(q;Pi,P2)(2Tr) 3 dq 2 , (174) 

giving the phase space factor in equation (171) as (see proof in Section 4.4), 

d®3(pi +P2-,P3,P5,P6 ) = d$ 2 {pi +P2]P4,P3)d$2{P4]P5,P6)(2n) 3 dpl . (175) 
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Substituting equations (170) and (175) into equation (171) gives 


da( 12 -> 356) = 


2 \ 2 |All2^34| 2 (2vr) 4 


4 P 1-^4^561 (jp\ - ml) 

X(M> 2 (j> 1 +P25P4,P3)rf < I > 2(P4;P5,P6)(27r) 3 (ip| 


= dcr(12 -»• 34) 


2 7774 
(2VT) 4 


dT(4 — > 56)- 


or 


dcr(12 — > 356) = d<r(12 — > 34) dr (4 — > 56) — - 

7 r 

Integration of da and dr gives 

r m.ri'A 

a(12 -> 356) = / dp^ <r(12 -> 34) 


(pi - m\) 2 



(176) 

2>2 d ‘Pl ■ 

(177) 

~ W.4)- 


56)/7t 

.2X2 

(178) 


De Wit [23] (p. Ill) and Pilkuhn [18] (p. 36) have taken into account the decay width of the 
intermediate A resonance by the substitution 


1 


1 


p 2 + m 2 p 2 + m 2 — imY 


Here, Y is distinguished from the T(4 — > 56) of equations (173) and (176) by 

r(4— >56) = partial width, 
r = total width, 


(179) 


(180) 

(181) 


where the total width includes all possible decay paths. However, the propagators used in the 
present work are defined with a different sign, and therefore the above substitution becomes 


1 


1 


p\ — m\ pi — ml — im^Y 

The square of a complex number is (x + iy ) 2 = x 2 + y 2 so that 

1 1 

{pi ~ ml ) 2 (pi - ml ) 2 + m 2 A Y 2 ’ 

The differential cross section in equation (177) becomes 


(182) 


(183) 


da (12 -> 356) = dcr(12 — 34) dT(4 -> 56) 


m4 


7T (pi - ml) 2 - (m A r) 2 


dpi , (184) 


and the total cross section in equation (178) becomes 


ct(12 -> 356) 


I dpi a (12 — > 34) 


777-4^4 — > 56) / 7T 
(pI - ml) 2 + 777. 2 r 2 


(185) 
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This can be re-written using [24] 


p{^ 2 ) = 


m 4 T / 7r 


(p 2 — m|) 2 + m\ r 2 


where p 2 = p\, to give 


a (12 -> 356) = / dp 2 <r(12 -> 34) p(^) 


2\ r(4^56) 


Relabeling p 2 = x, the integral above is simply 

f dx 1 . 

J (x — m 2 ) 2 + m 2 T 2 

See section 4.2.2 below for an alternate derivation of equations (184) and (185). 


(186) 


(187) 


(188) 


4.2.1 Integration bounds on p 2 

In equation (185), the upper integration bound on the value of p\ is imposed by y/s, the total 
center of momentum energy of the interaction. In the reaction 12 — ► 34, the square of the total 
center of momentum energy is given by 

s = (pi +P 2) 2 = (P3 + P4) 2 

= mjj + ml + 2£' 3 I?4 - 2p 3 • p 4 

= m\ + m\ + 2E 3 E 4 + 2\pf\ 2 , (189) 

where p 3 = — p 4 and |p/| = |p 3 | = | p 4 1 . Equation (189) shows that for a given s, m\ varies 
with |p/|, and has a maximum value when |pj| = 0, that is when all the mass and energy of 
particles 1 and 2 go into the mass of particles 3 and 4, with no kinetic energy remaining. In this 
case, E 3 = m 3 , E 4 = m 4 , and equation(189) becomes 


s = 

Solving for m 4 , and noting that p 2 = m 2 

(m 3 + m 4 ) 2 . 

(190) 

^4, max = 

Vs - m 3 , 

(191) 

3-max = 

pi = (Vs - m 3) 2 ■ 

(192) 


The lower integration bound on the value of p\ is imposed by the masses of the decay products 
of particle 4. In the reaction 12 — » 34 — ► 356, particle 4 decays by 4 — >■ 56 into two real particles. 
For particle 4 to decay into particles 5 and 6, the mass of particle 4 must be at least as great as 
the sum of the masses of particles 5 and 6. Thus, 

m 4 , min = rn 5 + m 6 . (193) 
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This result can also be demonstrated as follows. In the rest frame of particle 4, the square of 
the center of momentum energy is 

s = rnl = {p 5 +p e ) 2 

= m\ + ml + 2 E 5 Eq - 2p 5 • p 6 

= ml + ml + 2Ez,Eq + 2|ps| 2 , (194) 

where the last step follows from p6 = — P5. The minimum value of m 4 occurs when |ps| = 

|p6 1 = 0, in which case E§ = m§ and Eq = me and equation (194) reduces to 

ml = (m.5 + ?«6) 2 , (195) 

which is equivalent to equation (193). 


4.2.2 Alternate derivation of the cross section 

An alternate, and better, derivation is now given. Equations (184) and (185) may also be derived 
by starting with the Feynman diagram for the 3 - body final state (figure 26) and noting that 
particle 4 is represented by an internal line (see Feynman Rules). This is a virtual particle whose 
mass 7774 can vary. Since particle 4 decays, it also has a decay width T. The internal decaying 
particle is described by the Breit-Wigner propagator [8] (p. 101), [16] (p. 163), 


Propagator = -5 5 . (196) 

P A — 7772 — 7 7774 T 

The amplitude J\4 for the reaction 1 + 2— >4+5 + 6 is then determined by applying the Feynman 
rules to figure 26, yielding 


\M \ 2 = {g n NN oInaY 


{p\ — 777^) 2 + (777AT) 2 \t ~ 777 2 77 — 777; 


+ 


(197) 


Substituting the amplitude from equation (197) leads to equation (184), which may be expanded 
to 


where 


d< ™ 6= (iSi 1 " 12 - 3412 *) (siflSp 1 ^ 4 - 5612 " 5 ) 

777.4 dpi 

TT (pi - m 2 A ) 2 - (777aP) 2 


777. A 

r 

7772 
I P5 1 


A resonance mass, 

A resonance width, 
p\ / rn A i n general, 

5 + v /m 4 + m 4 + m| 


2(777 2 777g + 777 2 777 2 + 777.g777|) . 


(198) 
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The result for |ps| is taken from reference [10]. S 34 and £56 are the statistical factors associated 
with the sets of final state particles (3,4) and (5,6), respectively. For the reaction NN — > NA — > 
NNn, both factors are unity. The integrals over dt and (AI 5 are performed holding ra| constant. 
Since T is invariant (like rest mass, T is a rest width), the integral over di\ is most easily 
performed in the cm frame for particle 4 and yields a factor of 47 t. While the integral over 
dt is independent of p\, the integration bounds to and t n are functions of rn\{= p 2 ), so that 
integration over dt must precede integration over p 2 . 

4.2.3 Narrow width approximation 

Consider the approximation of a narrow width. A narrow width (T — ► 0) corresponds to a long 
lifetime (r — > 00 ). The narrow width approximation allows equation (186) to be simplified. The 
narrow width approximation gives [18] (p. 37) 

5(x — m 2 ) = 


so that equation (186) becomes 

<7(12 — > 356) = J dp 2 <r(12 — » 34) 5(p 2 — m 2 ) — — ^ ° U) , (200) 

to finally give the cross section in the narrow width approximation [18] (p. 37), 

ct(12 -> 356) = <7(12 -> 34) T<y4 y 56 ^ . (201) 

A further simplification is achieved if the partial width is nearly equal to the total width, i.e. 

T « T(4 56) . (202) 


lim p(x) 
mr^O 

lim 


mV / 7r 


mr ->0 (x — m 2 ) 2 + m 2 T 2 ’ 


(199) 


TtA 


Then 


cr(12 -»• 356) w a(12 -> 34) , 


(203) 


or 


a(NN -»• NN tt ) w a(NN NA) . (204) 

This expression states that if the A resonance decays primarily to a particular final state Nir, 
then the total cross section for 7 r production is nearly the same as the total cross section for A 
production. Since the calculation of the cross section of a 2 - body final state NA is simpler than 
that of a 3 - body final state NN n, equation (204) is a considerable simplification in calculating 
the total cross section for 7 r production. 


42 



4.3 dLips formalism 

In this subsection it will be shown how to develop useful formulas to handle Lorentz invariant 
phase space. Note that the above is discussed by [18] (pp. 16, 18) who uses a formalism involving 
the Lorentz invariant phase space dLips(pi, ...p n )- The result of Pilkuhn is easily shown to be 
identical to the result obtained here. We have the definition [18] (p. 16) 


dLips(pi . . . p n ) = J| 


d 3 Pi 


1 = 1 (2ir) 3 2Ei 


and 


dLips(s; pi . . . p n ) = (2ir) 4 5 4 (P - ^pi) dLips (pi . . .p n ) 

i 

n n ,3 

= (27r) 4 (M> n (P;pi . . .p n ) . 

The Pilkuhn recurrence formula is (with s d = q 2 = (p\ + P 2 ) 2 ) 

dLips(s-,pi,p 2 ,P 3 ) = ;z—dLips(s;p d ,p 3 )dLips(s d ;pi,p 2 )ds d 

Z7T 

= (2Tr) 4 d$> 3 (P-,pi,p 2 ,p 3 ) = ^-{2'K) 4 d$ 2 (P;p d ,p 3 ){2ir) 4 d<i>2{q;Pi,P2)dq 2 


2t r 


or 


d$> 3 (P-,pi,p 2 ,P3) = ( 2ir) 3 d$2(P;Pd,P3)d$2{q;Pi,P2)dq 2 


in agreement with equation (175). 


(205) 


(206) 


(207) 


(208) 


4.4 Proof of phase space formula 

We now prove the phase space formula in equation (175), 

d$ 3 (pi +P2]P3,P5,Pe ) = d$ 2 {pi +P2]P4,P3)d^2{P4-pP5,P6){^) 3 dpl . (209) 

The proof proceeds as follows [18] (p. 18). The various phase space factors are 

d 3 p 3 d 3 p 5 d 3 p§ 


d$> 3 (pi+p 2 ;P3,P5,P6) = < 5 4 (pi +P2 -P3 ~Ph ~Pe) 

d$ 2 (p\ +P2',P4,P3) = d 4 (pi+P2~P4-P3) 

d^2(p4;P5,Pe) = ^{p4 - P5 - Pa) 


(: 2it) 3 2E 3 ( 2ir) 3 2E b ( 2ir) 3 2E 6 ’ 
d 3 p4 d 3 p 3 


{2tt) 3 2E 4 ( 27 t ) 3 2£ 3 ’ 
d 3 p 5 d 3 p 6 


(: 2tt) 3 2E 5 ( 2ir) 3 2E 6 
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d 3 p refers to the volume element of a 3- vector and d 4 p refers to the volume element of a 4- vector. 
However, when we write p 2 or dp 2 , then the p refers to a 4- vector, not a 3- vector! This is why 
some authors [24] prefer to write dp 2 = dp 2 where p is the mass, given by p 2 = p 2 = E 2 — |p| 2 . 
Using the identity 


1 = d 4 p 4 5 4 (p 4 - p 5 - Po) , 


( 210 ) 


yields 


d<5>z{pi + P2' 1 P3,P5,P6) = $ 4 {pi +P2-P3-P5- Pa) ^n) 3 2E 3 

xd 4 p 4 5 4 (p 4 -p 5 -P6) 


d 3 p3 
k) 3 2 
d 3 p 5 d 3 p e 


(2tt) 3 2E 5 (2tt) 3 2E 6 


= 5 4 {pi +P 2 -P 3 -P 5 - Pe ) 

Now using the result 1 = 8(x — a)dx or 

1 = 5 (pi - s d )dp\ 

with [24] (equation 14), [18] (p. 15), 

Sd = (Pb+Paf 


d 3 p3 


(2vr) 3 2 E 3 


d 4 pi d$ 2 (PA]p3,P(>) ■ 


gives 


— (-®5 + E&Y — (P5 + P6) 2 t 


d 4 p 4 = d 4 p 4 5(pl - s d )dp\ 

d 3 P4 , 2 (n 13 d 3 p 4 , 2 
~dp 4 = (2ir) ———dp 4 


( 211 ) 


(212) 


2 E 4 


(2tt) 3 2U 4 


where we have used the result 


d 3 P c / 2 2\ A 

— = d(p — m )d p . 


(213) 


(214) 


Thus, the above phase space factor becomes 


d®3(pi +P2-,P3,P5,Pe) = S 4 (Pl+P2-P3-Pi) 


d 3 p3 d 3 p 4 


(2tt) 3 2E 3 (2tt) 3 2E 4 


y dp'i d$2(p4',P5,p(i) 


= dd> 2 (pi +P2;P4,P3)d^2(p 4 ;P5,P6)(^) dpi 


(215) 


which is the desired result. 

The equation (214) used in the preceding step is derived as follows. Differentiating p 2 = 
E 2 — p 2 gives 


d A p = dE d 3 p 


(216) 
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(217) 


Expanding S(p 2 — m 2 )crp gives 


Applying the formula 


S(p — m 2 )d A p = 5(E 2 — p — m z )dE d 6 p 


5(f(E)) = J2 J^) S i E ~ E i) 


(218) 


with f(E) = (E + \/p 2 + m 2 )(E — i/p 2 + m 2 ), and keeping only the positive E root, gives 


5(p 2 — m 2 )d A p = S(E — i/p 2 + m 2 )dE . 

2jC/ 

Integrating over dE cancels the delta function, leaving 

d 3 p <rr 2 2 \ j 4 

— = d(p — m )d p . 


(219) 


(220) 


5 Conclusions 

We have presented some applications of a scalar quantum field theory method by considering 
massive particle exchange, and we have calculated differential and total cross sections for elastic 
and inelastic processes. The use of Mandelstam variables makes the mathematical expressions 
much simpler. The cross section formulas are all obtained in closed form and are compact 
and simple. The cross sections display typical behavior of elastic and inelastic processes. The 
theory only includes the lowest order terms. At larger energies, higher order terms are expected 
to contribute. However, these complications may sometimes be circumvented by appropriate 
adjustment of parameters, such as the coupling constants. The theory developed herein will be 
useful in calculating cross sections for production of all types of hadrons, if the effects of spin 
can be neglected. These cross sections can then be included in space radiation transport codes. 
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